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1.0  INTRODUCTION  AND  THEORY 


1.1  BACKGROUND  AND  PREVIOUS  MEASUREMENTS 

The  production  of  condensed,  or  two-phase,  flow  by  means  of  a 
supersaturated,  supersonic,  or  hypersonic  expansion  has  been  observed 
by  many  investigators.  Because  of  an  obvious  desire  to  predict  the  on¬ 
set  of  condensation  in  aerodynamic  simulation  facilities,  Daum  and 
Gyarmathy  (Ref.  1)  have  attempted  to  correlate  (with  the  nozzle  source 
reservoir  parameters)  the  observed  pressure  and  temperature  of  both 
nitrogen  and  air  at  which  onset  occurs.  Of  a  more  fundamental  nature, 
Wegener  (Ref.  2)  has  studied  homogeneous  condensation  processes  in 
expansion  flows  with  the  goal  of  obtaining  an  understanding  of  the 
relationship  between  onset  and  subsequent  condensate  growth  and  molec¬ 
ular  parameters.  The  approaches  taken  by  the  authors  of  both  Refs.  1 
and  2  are  classical  in  the  sense  of  using  concepts  of  critical  radii  to 
describe  the  onset  and  growth  processes,  and  both  references  contain 
excellent  bibliographies  of  previous  work  using  this  approach.  There 
has  been  a  relatively  recent  revival  of  interest  in  this  particular  prob¬ 
lem  for  a  variety  of  reasons,  including  the  possibility  of  the  use  of 
massive  clusters  or  polymer  beams  for  fusion  applications,  as  well  as 
interest  in  the  possible  enhancement  of  scattered  intensity  of  radiation 
either  incident  upon,  or  emitted  try,  a  plume  containing  clustered  mole¬ 
cules.  As  a  result  of  this  interest,  numerous  studies  of  clustered  or 
condensed  molecular  beams  have  resulted  using  either  mass  spectro- 
metric  (Refs,  3  and  4)  detection  of  polymers  in  the  free-molecular, 
far-field  region  of  the  plume  or  electromagnetic  filters  sensitive  to  the 
charge -normalized  number  of  molecules  per  cluster  (Ref.  5),  Addition¬ 
ally,  application  has  been  made  of  both  the  electron  beam  fluorescence 
method  (Refs.  6  and  7)  and  laser  Rayleigh  scattering  method  (Refs,  7 
and  8)  for  diagnostics  of  two-phase  flow  fields,  thereby  obtaining  even 
more  detailed  information  concerning  the  spatial  location  of  the  onset 
of  condensation.  Finally,  Raman  scattering  (Ref.  9)  has  been  applied 
to  the  problem  for  the  measurement  of  monomer  gas  density  throughout 
the  flow  field,  including  the  region  of  the  onset  of  condensation.  In  view 
of  the  added  measurement  capabilities  available,  it  is  felt  that  the  addi¬ 
tional  experimental  parameters  to  be  obtained  by  their  use  will  provide 
more  stringent  tests  of  proposed  theories  of  homogeneous  condensation. 
Moreover,  more  detailed  and  extensive  measurements  of  the  condensation 
process  are  required  for  expansion  flows  in  general  and  sonic  orifice 
flows  in  particular  to  provide  information  concerning  the  scaling  laws  of 
the  phenomenon  with  respect  to  orifice  diameter,  D;  reservoir  pressure, 
PD;  temperature,  T0;  and  gaseous  specie.  It  should  be  noted  that  previous 
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results  of  mass  spectrometric  sampling  have  been  used  to  provide 
scaling  laws  of  sonic  orifice  flows  with  respect  to  D,  PQ,  and  T0.  How¬ 
ever,  the  very  definition  of  the  occurrence  of  condensation  is  not  clearly 
delineated;  e.  g. ,  Daum  and  Gyarmathy  used  the  Wilson  point  criterion, 
while  Golomb,  et  al.  (Ref.  3)  defined  the  onset  of  condensation,  as  ob¬ 
served  in  the  farfield  of  the  plume,  by  two  criteria:  (1)  the  occurrence 
of  mass  spectra  peaks  of  dimers  (which  essentially  coincides  with  the 
n-mer  peak)  and  (2)  a  break  in  the  flow  speed  curve  when  observed  as 
a  function  of  reservoir  pressure,  Wegener  (Ref.  2)  generally  has  em¬ 
ployed  the  departure  of  the  static  pressure  from  the  isentropic  prediction, 
and  Hagena  and  Obert  (Ref.  5)  prefer  to  scale  according  to  constant 
values  of  (N/Z).  Obviously,  the  definition  of  onset  requires  some  care 
and  must  be  chosen  to  reflect  the  ultimate  application  of  the  results,  be-- 
cause  all  expansion  flow  fields  which  possess  finite  supersaturation  will 
exhibit  some  degree  of  condensation,  even  if  only  in  the  form  of  dimer 
molecules. 

It  is  the  purpose  of  this  work  to  provide  new  information  concern¬ 
ing  nitrogen  (N2)  condensation  in  a  sonic  orifice  expansion  flow.  Specifi¬ 
cally,  it  is  desired  to  obtain  information  concerning  the  spatial  location 
of  onset  of  condensation,  estimates  of  the  mass  fraction  of  condensate, 
and  information  concerning  both  the  spatial  variation  of  the  condensate 
growth  and  the  dependence  of  these  preceding  items  on  reservoir  pres¬ 
sure.  Finally,  with  such  information  it  is  desired  to  compare  these  re¬ 
sults  with  the  far-field  sampling  techniques  discussed  in  Refs.  3  and  4. 
For  this  comparison,  it  is  of  special  interest  to  obtain  information  con¬ 
cerning  the  size  range  and  number  density  of  N2  clusters  and  to  attempt 
to  correlate  previous  mass  spectrometric  observations  with  actual 
changes  in  the  flow -field  properties  such  as  gas  temperature.  To 
accomplish  this  goal,  it  was  decided  to  employ  both  the  electron  beam 
fluorescence  technique  for  far-field,  low-density  gas  temperature  diag¬ 
nostics  and  laser  Rayleigh  scattering  for  observation  of  the  condensate 
growth.  These  flow  diagnostic  techniques  will  be  briefly  discussed  in 
the  two  following  sections. 


1.2  ELECTRON  BEAM  FLUORESCENCE  TECHNIQUE 

The  use  of  the  fluorescence  radiation  resulting  from  inelastic 
electron-molecule  collisions  for  the  measurement  of  gas  density  and 
temperature  in  a  hypersonic  flow  field  has  been  demonstrated  by  many 
investigators  (Refs.  10  through  12).  Briefly,  a  high  energy  beam  of 
electrons  (with  energy  in  the  range  of  10  to  50  kev)  produces  upon  impact 
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with  either  atomic  or  molecular  species  excited  electronic  states.  By 
radiative  decay,  these  states  produce  a  fluorescence,  the  intensity  of 
which  is  directly  proportional  to  gas  density,  assuming  no  collisional 
quenching  effects.  Furthermore,  in  the  event  of  electron -molecule 
collisions,  the  molecular  radiation  consists  of  vibrational -rotational 
band  structure  within  the  electronic  transition  system,  which,  upon 
spectral  separation,  yields  information  concerning  the  temperatures 
of  the  internal  energy  modes  or  the  energy  distribution  function.  As 
described  in  Ref.  11,  in  the  case  of  Ng  the  electronic  system  used  for 
such  temperature  measurements  is  the  N<g  First  Negative,  n|  (1-), 
and  the  (0,  0)  vibrational  band  is  spectrally  resolved  for  determination 
of  the  rotational  temperature,  Tp,  of  the  ground-state  N2  molecule. 
Clearly,  observation  of  rotational  temperatures  in  excess  of  those  to 
be  expected  for  the  uncondensed,  but  rotationally  relaxing,  N2  flow 
will  be  a  manifestation  of  the  heat  release  attributable  to  the  recombi¬ 
nation  process  of  condensation  and,  thereby,  a  direct  diagnostic  for 
the  onset  of  condensation.  It  should  be  noted  that  the  electron  beam- 
induced  fluorescence  spectra  of  the  (0,  0)  band  of  the  N  2  (1-)  system 

are  interpreted  using  dipole  selection  rules  for  the  electron-molecule 
ionization -excitation  process. 


1.3  RAYLEIGH  SCATTERING 


It  is  well  known  that  the  photon-molecule  interaction  can  result  in 
a  scattering  event  for  which  the  internal  energy  state  of  the  molecule 
remains  unchanged.  The  scattered  photon  is  said  to  have  experienced 
Rayleigh  scattering,  and  the  scattered  photon  energy  is  essentially  the 
same  as  that  of  the  incident  photon.  Assuming  a  laser  source  of 
intensity,  IQ,  and  wavelength,  X,  Fig.  1  shows  the  scattering  configu¬ 
ration  used  in  this  study.  It  is  to  be  noted  that  the  incident  beam  is 
linearly  polarized,  and  the  polarization  direction  of  the  detector  is 
aligned  parallel  to  that  of  the  incident  beam.  It  is  easily  shown 
(Ref.  13)  that  the  scattered  intensity,  Is,  from  a  collection  of  N  mole¬ 
cules,  each  of  polarizability  a,  is  given  by 


Ka  -  N  .  a2/A4 


(1) 


where  the  constant  factor  includes  all  optical  transmission  and 
sensitivity  factors.  As  is  known,  the  polarizability  a  for  a  spiriless 
diatomic  molecule  is  given  by  the  trace  of  the  polarizability  tensor  a 


Tr2  = -a  -  (l/3)(a,  +  2a2) 


(2) 
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+  2.  Scattering  Direction 


Figure  1.  Polarization  direction  diagram. 


where  a\t  02,  and  03  represent  the  components  of  polarizability  along 
the  body  axes  of  the  molecule.  To  appreciate  the  value  of  Rayleigh 
scattering  as  a  diagnostic  for  the  onset  and  growth  of  condensation,  one 
can  simplistically  imagine  the  scattering  molecule  to  be  a  dielectric 
sphere  of  radius,  a,  and  index  of  refraction,  m.  The  polarizability  a 
can  be  replaced  by  the  expression  (Ref  /  13) 


a  =  (m2  —  l)a^/(ni2  +  2) 


and  Eq.  (1)  becomes 


Ka  x  N  x  a  ,/A4 


(3) 


where 


Ka  =  [(m2-l)/(n,2*2)J2x  Ka 
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If  the  scattering  medium  is  characterized  by  a  polydisperse  size 
distribution,  Eq.  (3)  becomes 

V1.  =  <  Nfl6>  (4) 

where  <Na®  >  is  the  expectation  value  of  the  product.  One  notes  from 
Eq.  (4)  that  upon  onset  of  condensation  the  scatterer  grows  in  size  as 
clustering  proceeds,  and  the  scattering  signal  varies  with  the  sixth 
power  of  the  radius  of  the  polymer.  Therefore,  using  the  simplistic 
view  that  a  dimer  is  twice  as  large  as  a  monomer,  only  1.  5  percent 
dimers  in  the  flow  will  essentially  double  the  scattering  signal.  There¬ 
fore,  Rayleigh  scattering  is  a  very  sensitive  measurement  for  the 
existence  of  clusters  in  a  condensing  flow  field.  It  should  be  noted  that 
if  the  cluster  size  is  such  that  a/X  &  0.05,  the  above  description  is 
adequate  and  the  more  complicated  Mie  scattering  theory  is  not  re¬ 
quired  (Ref.  14).  Previous  studies  of  condensation  processes  in  an 
expansion  flow  indicate  scatterer  sizes  on  the  order  of  100  A  or  less, 
and  for  visible  laser  wavelengths  the  Rayleigh  equations  are  adequate. 


Returning  to  the  polarizability  description  of  the  scattering  process, 
for  a  collection  of  polymers,  the  largest  of  which  is  designated  as  the 
j-mer,  the  scattered  intensity  Eq.  (1)  can  be  written  as 


1  /I0  = 

S  0 


(Ka/\*) 


N 


A  1 


i 


(N/NjKa 


«Al>2 

J 


(5) 


For  an  expansion  flow  field  it  is  desirable  to  nondimens ionalize  by 

use  of  the  reservoir  density,  N0,  and  Eq.  (5)  becomes 


I„/L 

S  O 


[■ 


=  (Kfl  x  N'0  x  a\ A4XN  i'Ha)  |l  +  k  (X/N  jXa 


i=  2 


(6) 


For  small  values  of  the  mole  fraction  of  condensate 


N 


and 


Nj/N 


~  X 
i  =  Ai 


Vnt 


where  is  the  mole  fraction  of  the  i-mer. 
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Now 


=  (K„ 


N. 


x  a 


2,A4)(Nj/N( 


[l  +  |2  X|<«/«1>S] 


(7) 


For  a  sonic  expansion,  the  ratio  Ni/N0  is  given  by  the  Sherman- 
Ashkenas  equations  (Ref.  15)  if  the  expansion  is  isentropic.  Assum¬ 
ing  AS  =  0,  Ni/N0isa  function  only  of  the  axial  distance  x/D  and  the 
specific  heat  ratio,  7.  Therefore,  for  no  condensation  and  relaxation, 
the  axial  variation  of  IS/ID  is  calculable,  and  any  deviation  from  the 
predicted  variation  can  be  attributed  to  a  nonisentropic  rate  process, 
the  type  of  which  is  indicated  by  the  variation  of  the  deviation  with  P0. 
Assuming  the  deviation  from  the  (change  in  entropy)  AS  =  0  prediction 
to  be  due  to  onset  of  nonequilibrium  condensation,  it  is  seen  that  the 
difference  between  the  measured  and  predicted  values  yields  an  experi 
mental  measurement  of 

i  Xj(a./aj)2 

i=2  111 


If  one  considers  flow  fields  for  which  only  dimers  and  monomers 
exist  and  assumes  Ni  approximately  equals  the  total  number  density, 
the  equation  for  Is/lo  becomes 

I9/l0  =  (Ka  x  N„  x  a^XNj/N^l  +  X2  x  (a2/ai)2]  (8) 

so  that  values  of  Xg  can  be  obtained  from  measured  values  of  Is/I0 
only  if  an  estimate  of  ct2lot\  can  be  found.  Assuming  additivity  of  bond 
polarizabilities  and  recognizing  the  weakness  of  the  van  der  Waals 
binding  of  nitrogen  dimers  (Ref.  13),  it  seems  reasonable  to  assume 
that  02  is  approximately  equal  to  twice  ari.  Thus,  Eq.  (8)  becomes 

I. /I.  -  x  ",  x  O|/A*)(N,/N0Ml  +  4Xjl  (9) 

Consequently,  for  low  values  of  condensate  mass  fraction,  for  which 
Ni/N0  is  approximately  equal  to  its  AS  =  0  value,  the  scattering  signal 
excess  to  the  isentropic  prediction  is  an  approximate  measure  of  the 
dimer  mole  fraction,  X2. 
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2.0  EXPERIMENTAL  APPARATUS  AND  METHOD 


2.1  GAS  SOURCE  AND  VACUUM  CHAMBER 

The  flow  field  was  produced  by  a  tubular  sonic  orifice  source  of 
3-mm  orifice  diameter  and  0.  013-cm  thickness,  a  diameter-to- 
thickness  ratio  of  approximately  24.  The  source  was  constructed  of 
stainless  steel  tubing  of  approximately  1-cm  inside  diameter  and 
10 -cm  length.  Copper  cooling  coils  were  wound  around  the  source 
to  provide  temperature  variation,  by  means  of  flowing  either  liquid 
or  gaseous  N2  through  the  coils,  and  several  strips  of  Nichrome® 
were  wound  atop  the  cooling  coils  to  aid  thermal  control.  The  source 
was  designed  so  that  the  ratio  of  the  residence  time  within  the  source 
and  the  rotational  relaxation  time  was  several  orders  of  magnitude 
greater  than  unity  to  ensure  thermal  equilibration.  The  source  was 
provided  with  copper- constantan  thermocouples  and  a  calibrated  pres¬ 
sure  gage  for  T0  and  Pc  measurements,  respectively.  All  data  were 
obtained  using  prepurified  N2  with  a  stated  purity  of  99.  998  percent, 
and  particulate  filtering  was  accomplished  using  two  Millipore® 

0,  025-^m  filters  arranged  in  parallel. 

The  gas  source  was  mounted  on  a  three-dimensional  traversing 
mechanism,  the  description  of  which  is  given  in  Ref.  8,  as  is  the  de¬ 
scription  of  the  vacuum  chamber  itself,  The  traversing  mechanism 
has  an  accuracy  and  reproducibility  of  0.  013  cm  in  the  axial  direction. 
The  vacuum  chamber,  which  enclosed  the  motor-driven  traversing 
mechanism,  was  maintained  at  background  pressures  of  10“®  torr  or 
less  with  no  gas  flow  and  10  "5  to  10”4  torr  with  maximum  gas  flow, 
and  Fig.  2  shows  a  section  view  of  the  chamber.  The  low  background 
pressure,  Pg,  with  mass  flow  rates  exceeding  1  gm/sec  was  main¬ 
tained  by  liquid  N2  and  gaseous  helium  (He)  cryopumping  liners.  Back¬ 
ground  pressures  were  measured  using  calibrated  ionization  gages,  and 
the  ratio  Pq/Pb  was  or  greater  for  all  cases.  Fluctuations  in  PD 
and  T0  during  the  course  of  data  acquisition  were  each  on  the  order  of 
one  percent  or  less. 


2.2  ELECTRON  BEAM  APPARATUS 

Optical  and  electronic  components  used  in  the  electron  beam  mea¬ 
surements  are  shown  in  Fig.  3.  The  electron  source  was  Radio 
Corporation  of  America  (RCA)  model  VC2126  V4  electron  gun  mounted 
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Figure  2.  Vacuum  chamber  schematic. 
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Figure  3.  Optical  and  electronic  arrangement  for  electron  beam  measurements. 
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in  a  stainless  steel  housing  at  the  bottom  of  the  chamber,  and  the  beam 
was  directed  vertically  upwards.  The  beam  was  injected  into  the  cham¬ 
ber  and  flow  field  through  a  1.0-mm-diam  orifice  and  was  magnetically 
focused  to  provide  a  sufficiently  small  diameter  beam  source  at  the 
centerline  of  the  cell  and  flow  field.  The  chamber  centerline  was 
located  approximately  50  cm  from  the  injection  orifice.  A  15-cm 
diffusion  pump  with  a  liquid -N2- cooled  baffle  maintained  the  electron 
gun  pressure  at  less  than  2  x  10“^  torr.  The  electron  gun  power 
supplies  and  voltage  dividing  network  are  shown  in  Fig.  4.  All  mea¬ 
surements  used  a  potential  of  40  kv.  The  beam  collector  cup  was  simi¬ 
lar  to  that  described  and  shown  in  Ref.  16,  except  that  its  length  was 
34  cm.  The  distance  from  the  plume  axial  centerline  to  the  first  grid 
of  the  collector  cup  was  30  cm.  Vertical  alignment  of  the  gun  exit 
orifice,  collector  cup,  and  pitot  probe,  which  served  as  a  plume 
centerline  reference  point,  was  performed  with  a  plumb  bob. 

The  beam  fluorescence  was  focused  onto  the  20-pm-wide  entrance 
slit  ofoa  0.  75-m  Spex  scanning  spectrometer  with  a  1200 -groove/mm, 
3000 -A  blaze  grating  and  operated  with  equal  entrance  and  exit  slit 
widths.  The  25.  4-cm  focal  length  lens  was  apertured  to  a  diameter  of 
4.  5  cm  to  match  the  spectrometer  f-number.  Scan  rates  were  varied 
from  1  to  3  A/ min.  The  slower  rates  were  required  for  adequate  line 
peak  duration  times  in  low-density  flow  regions. 


Figure  4.  Schematic  of  electron  gun  power  supply. 
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During  the  course  of  the  measurements,  Tp  data  were  acquired 
with  both  the  spectrometer  slit  parallel  to  the  electron  beam  and,  by 
use  of  a  Dove  prism,  perpendicular  to  the  beam.  The  slit  heights  were 
1  to  3  mm  and  6  mm,  respectively,  and  the  respective  magnification 
factors  were  1.  6  and  1.  8,  To  the  accuracy  of  the  measurement,  no 
systematic  differences  in  the  results  for  Tr  were  noted  with  the  two 
configurations,  and  it  must  be  concluded  that  no  systematic  effects  of 
flow -field  gradients  are  present  in  these  results. 

A  Calibration  Standards  Corporation  power  supply  provided  1500  v 
to  the  EMI  6256S  photomultiplier  tube  and  1800  v  to  the  International 
Telephone  and  Telegraph  (ITT)  FW130  photomultiplier  tube.  The 
6256S  tube  has  an  "S'^Q)  spectral  response  and  is  designed  for  very 
low  light  levels  in  the  near  ultraviolet.  The  FW130  has  an  S-20  spectral 
response  and  is  also  particularly  suited  for  photon  counting  applications. 
Figures  5  and  6  show  signal-to-noise  ratios  as  a  function  of  discriminator 
voltage  for  a  6256B  tube  at  its  operating  temperature  of  -20°C  and  for  the 
FW130  tube  at  -30°C.  The  behavior  of  the  6256S  is  close  to  that  of  the 
6256B.  Discriminator  voltages  of  0,  6  v  for  the  6256S  and  0.  5  v  for  the 
FW130  were  selected  for  an  optimum  balance  between  count  rates  and 
sufficiently  high  signal -to -noise  ratios.  The  signal -to -noise  ratios  are 
greatly  improved  at  cold  temperatures. 

In  obtaining  a  temperature  from  a  scan,  no  spectral  sensitivity 
corrections  were  necessary  since  only  relative  peak  heights  are  re¬ 
quired  in  a  narrow  wavelength  interval. 

The  photon  counting  system  consisted  of  the  following  Ortec  instru¬ 
ments:  Model  454  timing  filter  amplifier,  Model  436  100-MHz  dis¬ 
criminator,  Model  441  ratemeter,  and  Model  715  dual  counter/ timer, 

A  50 -ohm  photomultiplier  output  load  resistor  and  maximum  ampli¬ 
fication  of  the  electronic  processing  system  were  used.  Two  methods 
of  data  recording  were  employed.  For  the  first  method  a  Digital 
Equipment  Corporation  PDP®-8F  computer  was  used  to  acquire  the 
output  data  from  the  counter /timer.  The  second  method  consisted  of 
a  Honeywell  Electronik  19  strip- chart  recorder  registration  of  the 
analog  output  of  the  ratemeter. 

The  PDP  data  acquisition  program  employed  included  provisions 
for  spectral  background  measurements,  measurement  of  R -branch 
peak  amplitudes  and  wavelengths,  real-time  calculations  of  rotational 
temperatures,  and  normalization  of  the  data  for  variation  of  the  elec¬ 
tron  beam  current. 
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Figure  5.  Signal-to-noise  ratio  versus  discriminator  volt-  Figure  6.  Signal-to-noise  ratio  versus  discriminator  volt¬ 
age  for  EMI  6256  photomultiplier  tube.  age  for  ITT  FW-J  30  photomultiplier  tube. 
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The  spectrometer  data  were  acquired  through  the  digital  input 
buffer  from  the  Ortec  photon  counting  equipment.  A  variable  sampling 
interval  was  available  and  chosen  for  compatibility  with  the  spectrom¬ 
eter  scanning  speed  and  signal  intensity.  The  line  wavelength  was  de¬ 
termined  from  the  time  measurement  which  was  read  into  the  computer 
along  with  the  spectrometer  data.  The  Ortec  B  counter  was  used  in  the 
0.  1-sec  clock  mode.  The  spectrometer  starting  wavelength  and  scan 
speed  were  manual  inputs  via  the  keyboard  such  that  the  wavelength  of 
the  individual  lines  could  be  identified. 


An  average  background  count  was  determined  try  count  accumu¬ 
lation  for  several  sampling  intervals  and  simply  dividing  by  the  number 
of  samples  taken.  A  background  wavelength  interval  was  specified  via 
the  keyboard.  Similarly,  the  wavelength  interval  for  a  temperature 
measurement  was  input  to  the  computer. 


The  line  amplitudes  were  determined  by  selecting  the  maximum 
value  that  existed  between  a  specified  number  of  monotonically  increas¬ 
ing  samples  and  a  similar  number  of  monotonically  decreasing  samples. 
The  number  of  samples  to  identify  a  positive  and  negative  slope  was 
variable  and  generally  chosen  as  three.  This  approach  was  taken  to 
isolate  false  peaks  attributable  to  noise.  Although  three  samples  for  a 
slope  identifier  occasionally  gave  false  peaks,  this  number  of  samples 
was  chosen  for  sensitivity  purposes.  Real  and  false  lines  were  identi¬ 
fiable  by  inspection  of  the  associated  wavelengths. 


For  the  purpose  of  obtaining  real-time  temperature  results  from 
the  spectral  data  ratios  of  individual  R -branch  rotational  line  intensi¬ 
ties  of  the  N*£(l-)(0,0)  band  were  used.  Theoretical  calculations  of 
the  electron  beam-induced  fluorescence  of  the  Nj(l-)(0, 0)  band  were 
performed  using  the  spectrometer  dispersion,  entrance  and  exit  slit 
width  values.  Ratios  of  the  rotational  line  intensities  of  the  K  =  1 
and  3  levels  and  3  and  7  levels  were  expressed  as  a  function  of  Tr 
over  the  Tr  ranges  of  10  to  70°K  and  50  to  ISO'K,  respectively.  The 
rotational  temperature  was  calculated  on-line  by  the  computer  using 
the  equations 
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where  R(l),  R(3),  and  R{7)  are  the  rotational  line  intensities  of  the 
transitions  from  the  K  =  1,  3,  and  7  rotational  levels,  respectively. 
Sufficient  computer  storage  was  not  available  to  use  multiple  line 
least  squares  curve  fits  for  temperature  determination.  However, 
the  above  curves  were  adequate  to  the  extent  that  temperature  errors 
less  than  2°K  were  obtained  when  compared  to  least  squares  calcula¬ 
tions. 

Background  corrected  peak  heights  of  the  emission  lines  were 
used  in  computer  program  TROT  {Ref.  16)  to  calculate  the  rotational 
temperature  from  the  Boltzmann  plot  as  described  in  the  previous 
reference  to  the  experimental  data. 

Using  analog  recording  at  each  position  in  the  flow  for  which  data 
were  acquired,  background  deflections  for  each  scale  position  of  the 
ratemeter  were  recorded  at  3918  A.  This  wavelength  location  was 
devoid  of  any  impurity  radiation  and  yet  only  4  A  from  the  3914  A 
(0,  0)  bandhead,  thereby  precluding  any  wavelength  dependent  effects 
on  the  background  measurement.  The  spectral  scan  was  begun  at 
approximately  the  bandQorigin  at  3910  A,  and  scanning  was  concluded 
at  approximately  3900  A. 

The  ratemeter  scale  factors  were  selected  to  obtain  maximum 
chart  deflections  for  each  line,  and  the  ratemeter  was  linear  over  the 
ranges  selected.  The  time  constant  associated  with  the  standard  devi¬ 
ation  selector  of  the  ratemeter  was  kept  equal  to  or  less  than  0.  5  sec. 
Because  the  time  constant  is  a  function  of  the  scale  factor,  the  two 
switches  were  changed  together.  As  in  the  PDP-8  analysis,  program 
TROT  was  used  to  obtain  the  temperature  from  the  resultant  line  peak 
heights.  A  computer  program  that  weighted  each  line  according  to  its 
intensity  and  calculated  a  rotational  temperature  was  developed  for 
this  latter  part  of  the  measurements.  Since  stronger  lines  had  larger 
signal -to -noise  ratios,  more  weight  was  given  to  them.  The  rotational 
temperature  from  a  five-line  weighted  average  was  generally  taken. 


2.3  RAYLEIGH  SCATTERING  APPARATUS 

The  optical  and  electronic  experimental  arrangement  for  the 
Rayleigh  scattering  measurements  is  shown  in  Fig.  7.  The  laser  source 
employed  was  a  Coherent  Radiation  Model  52-B  argon-ion  laser  operating 
at  a  wavelength  of  4880  A  and  a  power  of  0,  5  w.  A  1/2-wave  plate  was 
used  to  rotate  the  normal  vertical  polarization  of  the  laser  output  to 
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horizontal  polarization  (the  polarization  vector  is  into  the  plane  of 
Fig.  7).  A  quartz  lens  of  3.  81-cm  diameter  and  1-m  focal  length 
focused  the  laser  beam  through  a  quartz  chamber  window  and  a  series 
of  three  0.48-cm-diam  apertures  to  a  point  approximately  on  the  cham¬ 
ber  centerline.  The  final  aperture  was  located  30.  5  cm  from  the  laser 
beam  focal  region.  As  the  laser  beam  defocused,  it  entered  a  laser 
dump,  the  entrance  of  which  was  a  black  cardboard  tube  which  extended 
into  the  chamber  to  within  22.  5  cm  of  the  chamber  centerline.  The 
latter  stage  of  the  dump  was  a  black  anodized  aluminum  tube  capped 
with  a  black  anodized  conical  section. 


Photon  Counting 
Equipment 


Black  Flocking 


Water  -*■ 
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Figure  7.  Optical  and  electronic  arrangement  for  Rayleigh  scattering  measurements. 


The  focal  region  of  the  laser  beam  was  located  at  the  focal  point 
of  the  collection  optics  system  which  was  oriented  at  90  deg  with  re¬ 
spect  to  the  flow-field  axis  and  laser  beam.  The  scattering  volume 
was  viewed  through  a  glass  port  by  a  glass  lens  of  12.  7-cm-diam  and 
24-cm  focal  length.  The  collected  Rayleigh  scattered  radiation  was 
transmitted  through  a  small  piece  of  HN-22  Polaroid®  material 
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with  the  polarization  transmission  axis  usually  parallel  to  that  of  the  o 
laser  beam  and  then  through  a  4880-A  pass  interference  filter  of  10-A 
half-width  (see  Fig,  8  for  transmission  curve).  The  radiation  was 
focused  onto  a  0.  0508-cm-diam  aperture,  and  light  passing  through 
the  aperture  was  focused  by  a  small  lens  onto  the  photocathode  of  an 
ITT  FW130  photomultiplier  tube  mounted  in  a  thermoelectrically 
cooled  housing  and  cooled  to  -26°C. 


Figure  8.  Filter  transmission  curve. 

The  collection  optics  system  was  shielded  with  black  cardboard 
and  shrouded  with  layers  of  green  felt.  For  further  reduction  of 
stray  light,  a  collection  optics  viewing  dump,  consisting  of  a  piece  of 
black  cardboard  tubing,  was  located  at  the  bottom  of  the  test  chamber 
collinear  with  the  scattering. volume  and  the  center  of  the  collection 
lens.  All  optically  accessible  surfaces  inside  the  test  chamber. 
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including  the  sonic  orifice  source,  were  either  painted  with  a  flat  black 
coating  or  covered  with  a  black  flocking  material.  An  additional  25.4-cm 
length  of  black  cardboard  tubing  was  attached  to  the  interior  side  of  the 
collection  optics  viewport,  and  a  3.  8-cm-diam  aperture  was  mounted 
on  the  end  of  this  tube  to  restrict  the  field  of  view  of  the  collection 
optics. 

The  combination  of  input  laser  beam  optics  and  collection  optics 
provided  a  spatial  resolution  of  approximately  1  mm^. 

Photon  counting  was  used  to  process  the  photomultiplier  tube  output 
signal  which  was  proportional  to  the  scattered  light  intensity,  Ig.  The 
photomultiplier  output  was  coupled  through  a  50 load  resistor  to  an 
Ortec  photon  counting  system.  Counting  time  intervals  were  preset  for 
either  300,  100,  or  10  sec. 

The  discriminator  setting  was  such  to  maximize  the  signal-to-noise 
ratio,  and  the  overall  linearity  of  the  system  had  been  previously  veri¬ 
fied.  The  counting  rates  encountered  during  the  investigation  were 
sufficiently  low  relative  to  the  20-MHz  capability  of  the  system  that  no 
pulse  pile-up  effects  existed.  Finally,  linearity  of  the  signal  with  laser 
power  had  been  verified,  and  the  observed  fluctuations  in  incident  laser 
beam  intensity  were  negligibly  small. 

Alignment  of  the  laser  beam  focal  volume  with  the  focal  point  of  the 
collection  lens  was  necessary  before  actual  scattering  measurements 
began.  Initially,  the  small  diameter  pitot  probe  was  swung  into  a  posi¬ 
tion  that  placed  it  on  the  centerline  of  the  test  chamber.  A  small  incan¬ 
descent  light  with  frosted  bulb  was  placed  just  below  the  tip  of  the  pitot 
probe  with  the  flat,  top  surface  of  the  bulb  just  touching  the  underside 
of  the  probe  tip,  and  the  center  of  the  circular  bulb  face  was  located 
0.  318  cm  fromthe  probe  tip  and  collinear  with  the  probe.  Collection 
optics  were  adjusted  in  such  a  manner  that  the  image  of  the  lighted  bulb 
was  in  focus  and  centered  on  the  scattering  volume  defining  aperture 
(an  0.  0508-cm-diam  pinhole).  With  the  lightbulb  removed,  and  the 
laser  operated  at  low  power,  the  laser  beam  focal  volume  was  positioned 
0.  318  cm  from  the  probe  tip  and  on  the  axis  of  the  pitot  probe.  With  the 
laser  operating  at  high  power,  and  with  cigarette  smoke  blown  gently 
into  the  focal  region,  the  laser  beam  focal  volume  was  visually  aligned 
with  the  scattering  volume  defining  aperture  by  slight  angular  adjust¬ 
ment  of  the  collection  lens. 
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In  order  to  relate  the  measured  relative  scattering  intensity, 

(Is/I0),  to  values  of  <No,^>  and  N,  calibrations  had  to  be  performed. 

For  each  calibration,  the  test  chamber  was  evacuated  at  room  tempera¬ 
ture.  After  the  chamber  pressure  had  decreased  below  0.05  torr,  the 
laser  was  operated  at  0.  5  w  and  the  background  (or  zero  pressure) 
relative  intensity  was  observed. 

After  valving  off  the  chamber  vacuum  pumps,  nitrogen  was  slowly 
leaked  into  the  chamber  through  the  free -jet  gas  flow  system.  Rayleigh 
scattering  intensity  measurements  were  made  at  pressures  of  200  and 
then  400  torr,  as  indicated  by  a  calibrated  Wallace  and  Tiernan  gage. 

At  the  200-torr  point,  the  collection  lens  was  adjusted  slightly  to  maxi¬ 
mize  the  collected  radiation  passing  through  the  defining  aperture. 

The  calibration  signal  of  the  background  adjusted  Rayleigh  scattering 
was  linear  with  pressure  at  constant  temperature. 

Initially,  calibrations  were  performed  to  determine  the  constant 
Ka  in  Eq.  (5)  using  a  gas  sample  for  which  Ni  was  known  and  no  clusters 
existed.  Nitrogen  at  280  K  and  200  torr  was  used  for  the  N2  density 
reference  point.  However,  it  was  found  during  the  course  of  the  mea¬ 
surements  that  the  calibration  factor  determined  thusly  exhibited  varia¬ 
tions  from  day  to  day  as  a  result  of  small  disturbances  to  the  optical 
train  as  well  as  variations  of  position  during  the  evacuation  and  cooling 
of  the  vacuum  chamber.  Consequently,  it  was  decided  to  use  the  flow 
field  of  the  sonic  orifice  for  calibration  purposes.  To  do  so  requires 
that  the  axial  position  at  which  the  calibration,  or  normalization,  is 
performed  be  described  by  the  Sherman-Ashkenas  equation  (Ref.  15) 
for  N/N0  as  a  function  of  x/D,  that  the  flow  field  not  be  saturated  and 
that  no  rotational  relaxation  effects  be  present.  The  axial  position 
used  for  calibration  was  chosen  to  satisfy  these  criteria,  and  the  points 
so  selected  are  denoted  on  the  figures  presenting  these  results  as 
reference  points. 

Before  actual  flow -field  measurements  could  begin,  the  location  of 
the  free- jet  orifice  relative  to  the  observed  scattering  volume  had  to  be 
accurately  known.  The  axial  location  of  the  exit  plane  was  readily 
determined  by  axially  traversing  the  free  jet  in  0,  01-in,  (0.  0254-cm) 
steps  toward  the  laser  beam  path  until  the  side  of  the  free  jet  just 
touched  the  laser  beam  and  prevented  its  passage  across  the  test  cham¬ 
ber.  With  the  chamber  at  atmospheric  pressure,  the  axial  zero  position 
was  easily  observed.  The  vertical  and  radial  zero  positions  were  mea¬ 
sured  by  establishing  a  free.-jet  flow  field,  making  vertical  and  radial 
intensity  profile  measurements  and  determining  the  zero  positions  from 
the  flow -field  symmetry. 
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The  measurement  procedure  using  the  photon  counting  system  was 
as  follows: 

1.  The  free -jet  position  was  set  and  recorded, 

2.  With  no  gas  flow  and  the  chamber  evacuated  to  10"  ^ 
torr  or  less,  a  background  relative  intensity  measure¬ 
ment  was  made,  and 

3.  With  the  free -jet  flow  stabilized  at  the  desired  reser¬ 
voir  pressure  and  temperature.  Step  2  was  repeated, 
and  the  number  of  photon  counts  per  second  minus  the 
background  signal  gave  the  Rayleigh  scattered  relative 
intensity  in  counts  per  second. 


3.0  RESULTS  AND  DISCUSSION 


3.1  ROTATIONAL  TEMPERATURE  RESULTS 


Rotational  temperature  measurements  were  performed  using  the 
previously  described  electron  beam  fluorescence  technique.  The 
reservoir  pressure  was  varied  over  the  range  of  100  to  1000  torr  at 
300  K  reservoir  temperature,  and  both  axial  and  radial  spatial  surveys 
were  performed.  Figure  9  shows  the  axial  variation  of  the  100  torr, 
or  PqD  value  of  300  torr-mm,  results,  and  shown  for  comparison  are 
the  results  of  Marrone  (Ref.  12)  obtained  for  P0D  equal  to  250  torr-mm. 
The  isentropic  curve  from  Ref,  15  is  also  shown,  and  the  difference 
between  the  values  of  Tr/Tq  and  T/T0  has  been  shown  conclusively  by 
Marrone  to  be  attributable  to  rotational  relaxation  effects.  The  1500- 
and  3000 -torr-mm  axial  variations  are  shown  in  Fig.  10,  and  Figs.  11 
and  12  show  the  radial  profiles  of  T^/T0  foi  the  300-torr-mm  condition. 
Figure  11  shows  for  comparison  purposes  for  the  x/D  =  6  profile,  data 
reported  by  Marrone  for  x/D  =  6.3  and  F0D  =  250  torr-mm.  The  agree¬ 
ment  is  good,  although  there  is  an  increasing  systematic  difference  in 
results  as  the  radial  distance  increases.  This  difference  is  most  likely 
real  and  is  due  to  the  different  values  of  ratios  of  reservoir  pressure, 
PQ,  to  background  pressure,  Pq,  uspd  in  the  two  measurements.  The 
data  of  Ref.  12  were  acquired  for  the';!.range  200  $,P0/Pg  1500,  where¬ 

as  the  present  data  were  acquired  with  P0/Pb  ~  106.  Consequently, 
the  present  sonic  orifice  flow  field  most  likely  experienced  greater 
expansion  and,  therefore,  lower  temperatures  near  the  jet  ''boundary.  " 
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Figure  9.  Axial  profile  of  TR/T0  -  low  P0D. 


Figure  10.  Axial  profile  of  TR/TD  -  high  P„D. 
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Figure  11.  Radial  variation  of  N2  rotational  temperature. 


Figure  12.  Isometric  profile  of  Tr/Td  of  N2. 
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To  determine  the  value  of  P0  at  which  condensation  onset  occurred 
within  the  flow,  measurements  of  Tr  were  made  at  discrete  axial  posi¬ 
tions  over  a  wide  range  of  PD  values.  It  was  obvious  that  at  low  values 
of  PQ  rotational  relaxation  resulted  in  values  of  Tr/Tq  >  T/T0,  and  as 
P0  increased  the  effect  of  binary  collisional  energy  transfer  decreased 
this  difference  in  temperatures.  In  principle,  neglecting  all  other 
effects, 

Lim  <Tb/T„)  .  T/Tc 

r„  -*■  a© 


However,  at  some  value  of  PD  condensation  onset  occurs,  and  the  heat 
of  recombination  increases  the  value  of  Tr/Tq  once  more.  Therefore, 
one  should  observe  a  minimum  in  the  value  of  Tr/Td  when  observed  as 
a  function  of  P0. 

This  series  of  measurements  was  performed  at  x/D  =  15,  20,  and 
30.  Figure  13  shows  the  results.  The  asymptotic  limiting  values  of 
the  T/T0  curve  are  shown  at  P0D  =  10^  torr-mm.  Using  the  values  of 
Tr/T0  determined  for  the  uncondensed,  rotationally  relaxing  flow 
fields,  a  curve  was  extrapolated  to  the  high  pressure  limit,  which  is 
the  isentropic  prediction  of  Ref.  15.  Deviation  of  the  measured  value 
of  Tr  from  this  extrapolated  curve  is  considered  as  a  manifestation  of 
the  onset  of  condensation.  For  these  three  axial  positions  of  x/D  =  15, 
20,  and  30,  the  average  onset  pressure  was  found  to  be  for  the  3-mm- 
diam  orifice 


(P0)g  =  290  ±_  50  torr 

Data  were  acquired  for  the  axial  variation  of  Tr/Tg  at  T0  =  116  K 
and  PQ  =  100  torr  and  T0  =  133  K  and  P0  =  500  torr.  Figure  14  shows 
the  result.  *  Additionally,  radial  variations  at  x/D  =  6,  10,  and  20 
were  obtained  for  both  cases  for  the  purpose  of  investigating  any  over¬ 
all  plume  characteristic  changes  for  flow  fields  which  were  obviously 
condensed.  These  results  are  shown  in  Figs.  15  and  16. 


*The  data  of  Fig.  14  have  been  corrected  for  rotational  relaxation 
effects,  as  shown,  assuming  a  temperature -independent  rotational  re¬ 
laxation  time  and  source  Knudsen  number  scaling  of  the  300  K  data  of 
both  Ref.  12  and  the  present  work. 
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Figure  13.  Variation  of  TR/T0  with  reservoir  pressure. 


Figure  14.  Low  Te  axial  profiles  of  TR/T0. 
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Figure  15.  Radial  variation  of  N2  rotational  temperature,  PD  =  100  torr. 
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3.2  RAYLEIGH  SCATTERING  RESULTS 

Axial  variations  of  the  scattering  signal  Is/I0  were  obtained  for 
reservoir  pressures  500,  1000,  1500,  2000,  3000,  and  4000  torr. 
These  results  are  shown  in  Figs.  17  through  22,  and  Fig,  23  shows 
all  results  on  a  single  plot.  The  value  of  (Ni/N0)^g_g  as  predicted 
by  the  Sherman-Asbkenas  equation  (Ref.  15)  is  shown  in  each  figure. 
The  deviation  of  the  data  from  this  curve  is  evident  for  P0  >  1000 
torr,  and  indicates  the  axial  location  of  the  onset  of  condensation, 
{x/D)q.  The  agreement  between  the  data  and  the  predicted  values 
for  1  £x/D  <  (x/D)g  as  shown  in  Figs.  17  through  23  is  to  be  noted. 


Axial  Distance,  x/D 

Figure  17.  Rayleigh  scattering  for  P„  =  500  torr. 


V 

*y 


Axial  Distance,  x/D 


Figure  18.  Rayleigh  scattering  for  P0  =  1000  torr. 


i 


0  2  4  6  8  10  12  14  16  18  20  22  24 


Axial  Distance,  x/D 


Figure  19.  Rayleigh  scattering  for  PD  =  1500  torr. 
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Figure  21.  Rayleigh  scattering  for  P„  =  3000  torr. 
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Also  shown  in  Fig.  23,  by  dashed  line,  are  the  results  obtained  for 
the  isentropic  expansion  density  ratio  in  the  region  x/D  <  1  for  which  the 
Ashkenas -Sherman  equations  are  invalid.  Reference  15  presents  static 
pressure  and  Mach  number  data  for  the  transonic  region  of  the  N2  ex¬ 
pansion  from  a  sonic  orifice  source.  Those  data  show  that  the  Mach 
number  at  x/D  =  0  is  approximately  equal  to  0.  70.  Using  this  value  of 
M  and  isentropic  flow  relations,  one  finds  the  number  density  ratio  pre¬ 
dicted  at  the  orifice  to  be  0.  792,  and  this  value  is  consistent  with  the 
extrapolated  data  of  Fig,  23.  It  is  also  noted  from  Fig.  23  that  the 
isentropic  density  ratio  of  0.  634  corresponding  to  a  Mach  number  .of 
one  is  attained  at  x/D  s  0.  25,  and  it  is  concluded  that  the  axial  sonic 
point  of  the  orifice  is  located  at  x/D  =  0,  25,  Therefore,  continued  ex¬ 
pansion  from  x/D  -  0.  25  to  x/D  =  1,  0,'  according  to  Ref.  15,  is  required 
for  development  of  a  self-similar  flow  field  with  an  apparent  source  lo¬ 
cated  at  x/D  s  0.  40. 

Figure  24  shows  the  expansion  adiabats  plotted  in  the  P-T  plane 
along  with  the  N2  vapor  pressure  curve,  which  was  obtained  from 
Ref.  17. 


Using  the  results  shown  in  Figs.  17  through  22,  the  value  of  (x/D)g 
was  determined  at  which  condensation  onset  was  just  detectable.  These 
values  are  shown  in  Fig.  25.  Also  shown  in  Fig.  24  are  the  x/D  loca¬ 
tions  determined  from  Figs.  17  through  22,  for  which  the  value  of  Is/I0 
exceeded  by  approximately  10  percent  the  isentropic,  uncondensed  value. 
It  is  seen  from  Fig.  25  that  as  the  value  of  PD  is  decreased  the  onset 
location  moves  away  from  the  saturation  location  to  positions  farther 
downstream,  indicating  higher  values  of  the  super  saturation  ratio 
s  =  Pg/PQ. 


Figure  26  shows  the  variation  of  s  with  P0  for  the  two  previously 
mentioned  limits  of  condensation.  The  magnitude  of  the  estimated  de¬ 
grees  of  supercooling,  Tq  -  Ts,  can  be  obtained  from  Fig.  25  and  the 
isentropic  curve  of  Fig,  9.  For  the  3-mm-diam  orifice,  it  is  found 
that  T#  -  Ts  =  45  to  50  K  for  PQ  500  torr,  -  Ts  =  20  K  for 
1000  PQ  <  3000  torr,  and  Tg  -  Ts  ^0  K  for  P0  >  3000  torr.  Finally, 
it  is  noted  that  the  spatial  gradient  of  temperature  at  onset  for  the  higher 
values  of  P0  is  approximately  37  K/cnri  or,  in  terms  of  the  time  rate 
of  change,  dT/dt  ^2.6  x  106  K/sec. 
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Figure  24.  N2  P-T  diagram. 
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Figure  26.  Variation  of  supersaturation  ratios  with  reservoir  pressure. 


3.3  ANALYSIS  AND  DISCUSSION 

Considering  the  Rayleigh  scattering  results,  it  is  recalled  that  con¬ 
densation  onset  was  not  observed  for  PD  =  500  torr,  but  deviations  from 
the  predicted  axial  dependence  of  Is/Io  were  seen  at  1000  torr.  Using 
the  mass  spectrometric  results  of  Ref.  3  and  using  PqD  scaling  for  con¬ 
densation  onset,  as  indicated  by  the  dimer  peak,  it  is  estimated  that 
only  dimers  exist  in  the  flow  field  for  P0  =  1000  torr,  D  =  3mm.  There¬ 
fore,  using  Eq.  (9)  the  axial  variation,  of  the  dimer  mole  fraction  for 
this  case  was  determined.  Figure  18;-shows  the  result.  Because  of  the 
very  low  signal  levels,  little  weight  has  been  given  to  the  point  at 
x/D  =  9.  6,  Additional  justification  for  this  is  given  by  comparisons  of 
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IS/ID  for  PQ  =  1500  and  2000  torr  where  it  is  seen  that  the  1000-torr 
point  at  x/D  =  9.  6  is  abnormally  high,  A  similar  calculation  for  X2 
for  PQ  =  1500  torr  is  shown  in  Fig.  19,  but  it  is  cautioned  that  these 
results  represent  a  maximum  for  X2  since  it  is  most  likely,  using 
mass  spectrometric  results  {Ref.  3),  that  larger  polymers  coexist  in 
the  flow  field.  Figures  21  and  22  show  the  axial  variation  of  the  con¬ 
densate  scattering  signal,  and  it  is  believed  that  this  is  indicative  of 
massive  condensation.  From  Fig.  23  the  condensate  growth  and  its 
cessation  or  freezing  point  is  clearly  indicated  for  PQ  =  3000  and  4000 
torr. 

The  axial  variations  of  the  Tr/Td  data  were  analyzed  for  the  pur¬ 
pose  of  obtaining  estimates  of  the  condensate  mass  fraction.  As  has 
been  discussed  by  Sherman  (Ref.  18),  the  enthalpy  equation  for  a  one¬ 
dimensional  condensing  flow  field  can  be  written  as 

«p  T„  T  f  *  L  .  Cp  Tc  -  „“/2  (10) 


where 


g  «  mass  fraction  of  condensate 
L  =  latent  heat  of  condensate 
uc  =  flow  speed  in  condensing  gas 

Cp  —  specific  heat  of  gas 

It  has  been  assumed  that  g  is  sufficiently  small  that  the  specific 
heat  ratio,  yM  remains  constant.  One  finds 

■0  -  iacpT0  +  e  x  l  .  <ptc#  ■  (n) 


and 


u  =  {2f c  T  -  c„T]l^  without  condensation 

p  0  p  (12) 

and  uc/u  >  1  if  gL/cpT0  >  <TC/T0)  -  (T/T0). 

That  is,  augmentation  of  flow  speed  can  result  by  condensation.  As 
x/D  -*  *,  (Tc/T0)  -*■  (T/T0)  -*  0  and  uc  -►  <1  +  gL/h0)*/2,  where  hQ  is  the 
specific  energy  of  the  reservoir. 
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For  small  changes  in  u  the  ratio  of  Eqs.  (11)  and  (12)  gives 


(Tr/To)  -  O7T0)  *  g  x  L/h0-  2  Au/u  (13) 

where  Au  =  uc  -  u  and  Au  can  be  either  greater  or  less  than  zero,  and 
the  temperature  of  the  condensed  flow, is  assumed  to  be  given  by  the 
measured  value  of  Tr.  To  estimate  g  from  the  Tr  data,  it  is  assumed 
that  the  condensation  process  occurs  instantaneously  at  (x/D )g  and  the 
area  change  of  the  plume  across  the  assumed  infinitesimally  small  con¬ 
densation  zone  is  zero.  One  can  then  use  the  well-known,  constant 
area,  diabatic  flow  equations  (Ref.  19)  to  relate  parameters  across  the 
condensation  zone.  For  diabatic  flow  with  constant  area,  it  is  known 
that  heat  addition  drives  the  Mach  number,  M,  toward  unity  for  both 
subsonic  and  supersonic  flow.  Therefore,  for  M  >  1  at  onset,  the  flow 
speed,  according  to  this  model,  initially  decreases,  and,  following  a 
re- expansion  of  the  flow,  the  flow  speed  increases  and  can  exceed  the 
isentropic,  uncondensed  flow  speed.  To  calculate  Au/u  and  g,  the 
Tr/T0  data  were  extrapolated  toward  the  onset  location  (x/D)g,  as 
indicated  by  the  Rayleigh  scattering  results,  and  Fig.  27  shows  this 
extrapolation.  The  difference  between  the  extrapolated  and  isentropic 
curves  does  not  exceed  the  experimental  uncertainty  in  the  Tr/T0  data 
until  x/D  equals  approximately  5,  which  was  assumed  to  represent  the 
location  of  condensation  onset  as  manifested  by  the  increased  value  of 
Tr.  At  this  point,  it  was  assumed  that  onset  was  instantaneous  and  of 
constant  area.  To  calculate  the  values  of  Au/u  and  g,  the  following 
procedure  was  used. 

The  following  constant  area  diabatic  flow  equation  relating  the 
temperature  and  Mach  number,  T  and  M,  respectively,  before  conden¬ 
sation  to  their  values  following  condensation,  Tc  and  Mc,  was  used: 

VT  =  KM^KKM) 

where 

KM)  =  M2/(l  +  yM2)2 

and  both  T  and  M  are  assumed  to  be  given  by  the  isentropic  equations 
of  Ref.  15  at  (x/D)g.  The  degree  of  condensation  is  assumed  to  be 
sufficiently  small  that  7  is  unchanged. The  measured,  or  extrapolated, 
value  of  Tr  at  (x/D)$  is  assumed  to  accurately  approximate  Tc.  Equa¬ 
tion  (14)  was  used  to  graphically  determine  the  Mach  number  Mc,  and 
using  the  measured  value  of  Tr,  the  flow  speed  following  condensation. 


(14) 

(15) 
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uc,  was  determined.  Therefore  uc,  coupled  with  the  flow  speed  value 
u  preceding  onset,  as  determined  by  the  isentropic  values  of  M  and  T, 
determines  the  value  of  Au/u  and,  thereby,  g,  using  Eq.  (13).  Addi¬ 
tional  diabatic  flow  equations  as  given  in  Ref.  19  were  used  to  graphi¬ 
cally  determine  the  values  of  gas  pressure  and  density  following  con¬ 
densation.  The  following  set  of  values  of  these  aforementioned  param 
eters  were  calculated  according  to  this  procedure: 


Before 

Condensation 


After  Final -to -Initial 

Condensation  Value  Ratio 


M 

T,  K 
N/N0 

u,  cm/ sec 
P,  torr 


6.  30 
33.  6 

4.30  x  10-3 
1.  182  x  104 
0.465 


6.  0 

36.  9 

4.  29  x  10'3 
1.  179  x  104 
0.  510 


0.  952 
1.  099 
0.  997 
0.  998 
1.  096 


The  mass  fraction  of  condensate  was  found  using  these  values  to  be 
approximately  1.4  x  10"2.  The  value  of  X2  from  the  results  of  Fig.  18 
is  approximately  1.  7  x  10" 2  or  a  mass  fraction  of  approximately 
3  x  10"2.  The  agreement  between  these  two  independent  determinations 
of  g  is  within  approximately  a  factor  of  two  which  is  perhaps  the  best  to 
be  expected  considering  the  simple  model  and  extrapolation  employed. 


Since 

Lim  uc  =  (1  * 
x/D  ■+  00 


and  such  velocity  augmentation  has  been  reported  by  Bailey  (Ref.  20), 
it  is  worth  noting  that  for  g  ^  1.  4  x  10”2,  the  magnitude  of  Au  of  N2, 
which  may  exhibit  rotational  relaxation  effects,  is  marginal  for  accu¬ 
rate  detection  by  time -of -flight  systems. 

It  has  been  noted  that  the  F2D  scaled  results  of  Ref.  3  would  indi¬ 
cate  a  predominance  of  dimers,  and  this  was  the  primary  basis  for 
proceeding  with  the  previous  analysis.  The  data  of  Yealland  et  al. 

(Ref.  4)  also  support  this  assumption.  Their  mass  spectrometric  re¬ 
sults  for  condensing  flow  fields  of  N2  were  fitted  to  a  functional  form 
relating  the  percentage  of  dimers  in  the  far  field  to  the  reservoir  con¬ 
ditions.  Using  P^D  scaling  and  the  equation  of  Ref.  4,  it  was  found 
that  Yealland  et  al.  predict  0.  08,  0.  56,  and  1,  92  percent  dimers  for 
PQ  =  500,  1000,  and  1500  torr,  respectively,  for  a  3-mm  orifice.  The 
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uncertainty  of  the  results  of  Ref,  4  is  stated  to  be  approximately  a 
factor  of  2,  and  these  results  are  in  order  of  magnitude  agreement 
with  the  present  results. 


Figure  27.  Axial  extrapolation  of  Tr/T0  for  diabatic  flow  analysis. 


A  similar  analysis  of  the  electron  beam  data  for  the  cases  of 
T0  =  116  K  and  133  K  was  performed.  The  T^/T0  data  were  extrapo¬ 
lated  to  the  saturation  point  <x/D)s  which  was  assumed  to  be  the  onset 
location  for  these  temperatures.  The  mass  fraction  of  condensate  was 
found  to  be  on  the  order  of  4,  0  x  10“^  and  1.  0  x  10 "2  for  T0  =  116  K 
and  133  K,  respectively.  However,  the  data  scatter  and  the  slight 
curvature  of  the  T0  =  133  K  data  extrapolation  give  greater  reliability 
to  the  lower  temperature  expansion. 


Concerning  the  possibility  of  gros's  changes  in  jet  characteristics 
of  condensing  flow  fields,  a  comparison  of  the  radial  profiles  of  Fig.  11 
and  Figs.  15  and  16  show  that  essentially  identical  radial  profiles,  non- 
dimensionalized  by  the  axial  value,  are  obtained.  However,  changes 
in  the  value  of  the  absolute  temperature  were  dramatic  for  flow  fields 

t, 
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with  condensate  mass  fraction  exceeding  approximately  one  percent, 
and  changes  in  T  were  detected  even  for  flow  fields  containing  conden¬ 
sate  in  the  form  of  only  dimers.  Insofar  as  previous  mass  spectro- 
metric  results  are  concerned,  these  results  aid  somewhat  in  quantifying 
these  previously  qualitative  results. 

It  should  be  noted  that  the  degree  of  supersaturation  predicted  by 
Fig.  1  of  Ref.  1  for  N2  flows  is  greater  than  that  actually  found  in  this 
work,  indicating  both  an  increase  in  sensitivity  for  the  detection  of 
condensation  onset  as  well  as  the  use  of  a  different  criterion. 


4.0  CONCLUSIONS 


The  use  of  both  the  electron  beam  and  Rayleigh  scattering  tech¬ 
niques  has  provided  independent  estimates  of  not  only  the  spatial  loca¬ 
tion  of  condensation  onset  hut  also  estimates  of  the  mass  fraction  of 
the  condensate  for  a  sonic  orifice  expansion  of  h'2-  The  mass  fraction 
of  condensate  results  are  in,  at  best,  order -of -magnitude  agreement 
with  estimates  of  previous  mass  spectrometric  results,  but  it  should 
be  noted  that  orifice  diameter  and  pressure  scaling  laws  were  used  in 
regimes  not  observed  try  the  far-field  mass  spectrometer  investigations. 
Additional  measurements  of  sonic  orifice  expansions  produced  by  large 
orifices,  i.  e. ,  millimeter  size  and  greater,  are  required  to  verify  the 
far-field  scaling  laws  and  to  observe  and  correlate  the  axial  growth 
variation  of  condensate.  More  detailed  information  concerning  the  gas 
temperature  is  required  near  and  at  condensation  onset,  and  measure¬ 
ments  in  progress  are  utilizing  rotational  Raman  scattering  for  this 
purpose. 

It  should  be  noted  that  the  degree  of  supersaturation  attainable  for 
the  high-purity  N2  flow  fields  studied  was  found  to  be  significantly  less 
than  that  predicted  by  Ref.  1,  and  it  is  suggested  that  the  recommenda¬ 
tions  of  Ref.  1  be  used  with  caution  in  the  design  of  uncondensed  N2 
hypersonic  flow  facilities. 
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NOMENCLATURE 


a  Particle  radius 

Cp  Specific  heat  at  constant  pressure 

D  Orifice  diameter 

g  Condensate  mass  fraction 

IQ  Incident  laser  beam  intensity 

Ig  Scattered  laser  beam  intensity 
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Ka  Optical  calibration  constant  for  Eq.  (3) 

Ka  Optical  calibration  constant  for  Eq.  (1) 

L  Latent  heat  of  sublimation 

M  Mach  number 

m  Index  of  refraction 

N  Number  density 

Number  density  of  i-mer 
N/Z  Number -to- charge  ratio  for  cluster 

P  Pressure 

Pj3  Background  pressure 

S  Entropy 

s  Supers aturation  ratio 

T  Temperature 

Tp  Rotational  temperature 

u  Flow  speed 

Mole  fraction  of  i-mer 
x  Axial  distance 

y  Radial  distance 

a,  a  Scalar  and  tensor  molecular  polarizability,  respectively 

Polarizability  of  i-mer 
AS  Change  in  entropy 

Au  Change  in  u  caused  by  condensation 

y  Specific  heat  ratio 

X  Wavelength 

SUBSCRIPTS 

v 

c  Variable  evaluated  for  condensing  flow 

o  Reservoir  value 
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s  Saturation,  value 

T  Total . 

6  Condensation  onset  value 

CGS  units  used  unless  stated  otherwise. 
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